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Abstract. A numerical transfer matrix calculation is presented for two fully- 
frustrated three-state Potts models on the square lattice: the Potts piled-up-domino 
model and the Potts zig-zag model. The ground state entropies and phase diagrams 
are found. The Potts piled-up-domino model displays a finite-temperature transition 
when the frustration effects are maximal, and displays reentrant behaviour, in contrast 
to the Ising model equivalent. The ground-state entropy per spin is larger than in the 
Ising equivalent. The Potts zig-zag model displays the same qualitative behaviour as 
its Ising counterpart, and the ground-state entropy per spin is the same in the Potts 
and Ising cases. 
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1. Introduction 

Spin models are said to be frustrated when they contain competing interactions which 
prevent the simuhaneous minimisation of all the lattice bonds in the ground state 
The system is additionally described as fully frustrated if each plaquette of the lattice 
is itself frustrated 10]. The nature of a phase transition depends partly on the symmetry 
group of the candidate low-temperature phases. Frustration effects are then of most 
interest when they lead to infinitely degenerate or "nearly-degenerate" ground states. 
These will in general satisfy a larger symmetry group than the individual spins, possibly 
giving rise to exotic critical phenomena One such phenomenon is re-entrance 0. This 
occurs when a disordered paramagnetic phase comes to lie at a lower temperature than 
an ordered phase; as the temperature is increased there is a transition from disorder to 
order. One possible explanation for this phenomenon is the so-called "order by disorder" 
mechanism 1^. 

In this article we present results for two fully-frustrated three-state Potts models: 
the Potts Piled-Up-Domino Model (from here on referred to simply as the piled-up- 
domino model) and the Potts Zig-Zag Model (or more concisely the zig-zag model). 
These are natural extensions of the Ising-model equivalents studied previously by Andre 
et a/.0 and are defined as follows: 

The g-state Potts model is defined as a set of variables {cXi} associated with the 
sites {«} of a lattice]]^. Each (Xj takes one of q distinct values. The Hamiltonian of the 
model is given by: 

where the sum runs over all pairs of sites and Jij are the interaction strengths between 
a given pair of spins. For the models of interest in this article, the interaction strengths 
are given by 

{0 for non-nearest-neighbour spins, 
Ji for n.n. spins along the solid lines, (2) 
J2 for n.n. spins along the dashed lines, 

as shown in figures |l] (A) and (B) for the Piled-Up-Domino and Zig-Zag models 
respectively. For convenience we define a = J2/J1. Ji will be taken as positive 
throughout. In this article we shall concentrate on the g = 3 Potts models. 

In the Ising versions of these models, the frustration effects suppress the critical 
temperature, which becomes zero when a = — We show that, while this remains 
the case for the three-state Potts zig-zag model, the three-state Potts piled-up-domino 
model has a finite critical temperature when a = — 1 and displays reentrant behaviour. 

The article is divided as follows: In the next section we present briefly the transfer 
matrix method used in our study, and in section |] we apply this method to the 
calculation of the residual ground-state entropy for the piled-up-domino model. We 
show that the residual ground-state entropy for the zig-zag model is exactly related to 
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the equivalent Ising-model case. In section ^ we present the finite-temperature phase 
diagrams for the two cases and section |^ will be devoted to discussion. 

2. The Transfer Matrix Method 

The results presented in this article are obtained using a (numerical) transfer matrix 
calculation, giving the partition function for strips of finite width L but infinite length, 
—>■ oo. The results are therefore numerically exact for a given lattice width L. 
The partition function for the model may be written in the form 

^ = Enm)^(C.;C.+i), (3) 

where V{Cx) is the contribution to the Boltzmann weight of a given column spin- 
state Cx due to interactions within a column and H{Cx] C^+i) is the contribution due 
to interactions between adjacent columns, and depends on the spin states of the two 
columns and C^+i. If additionally periodic boundary conditions are taken in the x- 
direction, the partition function for the system may be written as a trace over a product 
of transfer matrices: 

Z = Tr{T^}, (4) 

where T = V'^/'^HV^/^ and the sum over spin states is taken care of by the matrix 
multiplications. Denoting the eigenvalues of T by Aj, 

2 = EAf- (5) 

i 

As becomes large the partition function is dominated by the largest eigenvalue. 
The (dimensionless) free energy per spin, / is given by 

/ = lim \ogiZ), 

= ^log(Ao), (6) 
and the correlation length is given by 

where Aq is the largest eigenvalue and Ai is the second largest eigenvalue (in absolute 
value) 0. 

The transfer matrix T is indexed by and Cx+i both of which correspond to one of 
3^ distinct column spin states. The capacity required to store the 3^^ matrix elements 
on a computer soon becomes prohibitive. In order to extend the range of values of L 
which may be treated numerically, the transfer matrix T is decomposed into a product 
of sparse matrices as follows: 

One may view the transfer matrix as an operator which adds an additional column 
to the lattice. Each element of the matrix is indexed by the spin configuration on the 
last column of the lattice [L spins) and the spin configuration of the column being 
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added. If we define a dangling bond as the bond wliich has only one spin attached 
to it, we see that the transfer matrix replaces the L existing dangling bonds with L 
new ones. The entry configurations and exit configurations are given by the values of 
the spins attached to these dangling bonds. If now, instead of adding a whole column 
at once, one spin is added at a time, there are still L dangling bonds beforehand and 
afterwards (see figure This operation may be written as a matrix indexed by the 
configuration of the L spins connected to a dangling bond before the addition of the 
new spin, and the configuration of the L spins connected to a dangling bond after the 
spin addition. Since L — 1 spins are common to the input and output configurations, 
the matrix elements will be non-zero only if these spins have the same values in the 
input and output configurations. The corresponding matrix is thus sparse. In this way 
T may be decomposed: T = 7i7^^~^73 where the matrices 7i, 72, ^3 correspond to 
one-spin additions. If the boundary conditions and in-column interactions are taken 
care of when the first and last spins are added, then the intermediary spin additions 
may all be performed by the same matrix. The combination of the % being sparse and 
the small number of distinct matrices required to build up T vastly reduces the storage 
requirements and enables calculation for lattice widths up to L = 12. 

An alternative method for writing the transfer matrix involves the Kasteleyn- 
Fortuin mapping [|10| to map the Potts model to a lattice animal problem [|TT|. This 



mapping requires periodic boundary conditions to enable a sparse-matrix reduction |[TT 
For q = 3 the transfer matrices in the two representations require an equivalent amount 
of storage but in a spin representation it is easier to implement the different boundary 
conditions used in this article. 

The eigenvalues used were calculated using the Lanczos method, using the 
ARPACK package. 

3. Ground-State Properties 

In this section we calculate the ground-state entropy for the two models. In both 
models, when a > the model is ferromagnetic. While they are anisotropic, there is 
no competition in the interactions. The system minimises its energy by choosing the 
same state for all the spins giving q distinct ground states0, [1^. The entropy per spin 
is zero. 

A lattice bond is "satisfied" when the interaction energy along it is minimal, and 
"frustrated" otherwise. It may be seen that for a < it is no longer possible to 
simultaneously satisfy all four bonds of a plaquette. The model is then said to be 
frustrated. However, if a > — 1 there is no ambiguity, in the ground state, over which 
bond to frustrate. The frustration does not change the ground state configuration, and 
the entropy per spin remains zero. 

For a = —1 the energy of a plaquette is minimal regardless of which bond is 
frustrated as long as one and only one bond is frustrated. There are an infinite number 
of frustrated bond configurations satisfying the ground state condition and the entropy 
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per spin in this case is finite. 

For a < — 1 the ground state degeneracy is still infinite. While it is clear that the 
system will satisfy the antiferromagnetic bond in order to minimise the energy, there is 
still a choice over which of three ferromagnetic bonds to frustrate. The entropy per spin 
here is also finite. 

3.1. Piled-up- domino model 

The ground-state entropy for the piled-up-domino model was calculated using transfer 
matrices. In the ground state it is necessary to add two spins at a time in order to 
ensure that only one bond per plaquette is frustrated. Since the energy is constant and 
minimal, the entropy per spin is simply the number of ground states per spin, leading 
to: 

s=^log(Ao). (8) 

The entropy was calculated with transfer direction parallel to the dashed lines as shown 
in figure |l|. Both free and periodic boundary conditions were considered. The results 
are shown in table |1]. It may be seen that the entropy per spin is substantially larger 
than for the Ising-piled-up-domino model. 

The extrapolations to L — oo throughout the article are done using the Bulirch 
and Stoer algorithm |ll3[ . 

3.2. The zig-zag model 

The ground-state entropy per spin for the zig-zag model is the same in the Potts and 
Ising cases. This can be illustrated as follows: Imagine that, from some suitable 
boundary configuration of spins, a ground-state configuration of spins is built up by 
adding alternately plaquettes of type A and B, shown in figure |^, in rows, starting from 
top left and ending bottom right, one plaquette at a time. The spins marked by a cross 
are shared with plaquettes already added to the lattice, and are thus already fixed. In 
the ground state only one bond per plaquette is frustrated. Whichever of the bonds 
is frustrated in a plaquette of type A, the "free" spin is connected by a ferromagnetic 
path to a predetermined spin. At first sight this is not true for plaquettes of type B, 
however the frustrated bond is shared between two plaquettes, and it can be seen that 
the "free" spin is connected by a ferromagnetic path to a predetermined spin either 
directly, or via the neighbouring plaquette (of type A) see figure |^. The entropy per 
spin in the ground state is thus determined uniquely by the configurational entropy 
of the frustrated bonds, with no additional (local) degeneracy due to the number of 
Potts states. This was not so in the piled-up-domino case where it is possible to create 
free domains of spins. The residual ground-state entropy per spin is thus the same as 
reported for the Ising model[^, i.e.: 



s = — - / d6 d0 log 1 + 4 cos 6cos(j) + A cos^ 6 

4:TT^ JO Jo ^ 
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= 0.1615- ■■ for a< -1 

s = — / de d0 log 4(cos^ e + sin^ 4 

4:71^ Jo Jo L 

G 

= - = 0.2916--- for a = -1 

7r 

The difference between the ground-states of the Ising and Potts versions of the 
zig-zag models is an overall degeneracy relating to the symmetry of the spin {Z2 and 
Z3). This was not the case for the piled-up-domino model, where the entropy per spin 
for the Potts case was significantly larger than the for Ising case for both a = —1 and 
a < — 1. This difference indicates that the difference in symmetry group for the ground- 
states in the two cases is not simply due to the difference in the symmetry groups of 
the spins. Based on these results, it is reasonable to expect that the phase diagram of 
the Potts zig-zag model should be qualitatively similar to the Ising case, but that the 
finite-temperature behaviour of the piled-up-domino model may well be qualitatively 
different. This is indeed what will be observed in the next section. 



4. The Finite- Temperature Phase diagrams 

In this section we present the phase diagrams for the piled-up-domino and zig-zag 



models. The critical lines may be identified using phenomenological renormalisation fl^ . 
At a critical point the correlation lengths for two strip widths, measured as a fraction of 
the strip width, must be the same. This reflects the scale invariance of a critical system. 
In practice a finite width system is always off-critical, however 

L ~ U ^> 
gives a finite-size estimate of the critical temperature, T*, which will tend to the true 
critical temperature as the strip widths tend to infinity. Assuming that the system is 
sufficiently close to the true critical point, the correlation length behaves (to leading 
order) as 

iL = A\T* - Tc\~\ (10) 

Since it is the finite width of the system which prevents the correlation length from 
diverging, at T* we have L oc If these two scaling laws are admitted, then a finite- 
size estimate of the correlation length exponent may be calculated from the equation: 

1 log f ^ 



^L,L' log 



-1. (11) 



4.1. The zig-zag model 



The phase diagram found using phenomenological renormalisation for the zig-zag model 
with periodic boundary conditions is shown in figure ^. The finite-size effects are 
very small. The critical temperature is known exactly for two points on the phase 
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diagram ||T2|, |T^: a = 1, corresponding to the pure ferromagnetic Potts model on the 



square lattice, and a = 0, corresponding to the ferromagnetic Potts model on the 
hexagonal lattice. Our numerical values calculated at these two points are extrapolated 
and compared with the exact values in table 0. 

Qualitatively the phase diagram is the same as for the Ising model: There is a phase 
transition between a low temperature ferro-magnetic phase and a high temperature 
phase for a > —1. The critical temperature becomes zero at a = — 1 and there is no 
evidence of a phase transition for a < —1. 

The finite size estimates for the exponent u are shown in figure |[ For a pure three- 
state Potts model (a = 1 corresponding to the square lattice, and a = corresponding 
to the hexagonal lattice), the value of z/ = 5/6 is known exactly fl^. For a > the 
estimates for u converge to the exact value. Close to a = — 1 the estimates for i> seem 
to take on a different value, but the range of values of a over which this variation is 
observed diminishes as the system size increases. There is no evidence of a multicritical 
point in the phase diagram between a = 1 and a = — 1, and so it would seem reasonable 
to assume that the whole line up to a = —1 remains in the same universality class as 
the ferromagnetic three-state Potts model. 

4-2. The piled-up- domino model 

The phase diagrams for the piled-up-domino model with periodic and free boundary 
conditions found using phenomenological renormalisation are shown in figure M. 



The critical temperature is known exactly for a = 1 [|T2[, the pure three-state Potts 
model on the square lattice, and at this point the critical temperature coincides trivially 
with that found in the zig-zag model. It may also be easily determined for a = 0. When 
a = the spins on the dashed hues (see figure ^ (A)) may be traced over. The resulting 
model is a ferromagnetic Potts model on a square lattice with different interactions in 
the a; and y directions: = Ji and 

exp(2/j Ji) + 2 
^"P^^^^) = 2exp(M) + l- 
The critical temperature is known to satisfy 0]: 

(exp (/3 J,) - 1) (exp (pjy) - 1) = 3, 

which leads to = l/log(4), measured in units of Ji/k. In table ^ we compare 
the finite-size critical-temperature estimates, extrapolated using the Bulirch and Stoer 
algorithm [|l^, with the exact value calculated at a = 0. The agreement is good. 



Before trying to interpret the phase diagrams shown in figure ^ we review the 
main features of the Ising piled-up-domino model, reported in 0. The phase diagram, 
shown in figure |^, is made up of two transition lines in the Ising model universality 
class. The transition line for a > — 1 corresponds to the standard Ising transition 
between a ferromagnetic low-temperature phase and a paramagnetic high-temperature 
phase. For a < —1 the ground state is highly degenerate. The low temperature phase is 
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characterized as follows: the spins on the antiferromagnetic lines order ferromagnetically 
in the y direction and antiferromagnetically in the x direction, where x and y directions 
are taken as in figure |l|. The spins on the horizontal ferromagnetic lines are partially 
disordered. The transition is between this partially ordered antiferromagnetic phase and 
the high-temperature paramagnetic phase. The two transition lines meet with infinite 
slope at a = —1, where the critical temperature is zero. 

These results contrast with the phase diagrams shown in figure |^ for the three-state 
Potts models. The phase diagrams shown are calculated using finite width strips of width 
L in the direction perpendicular to the antiferromagnetic lines and infinite in the parallel 
direction. In figure ^ (A) periodic boundary conditions are taken and in figure ^ (B) 
free boundary conditions are taken. Periodic boundary conditions reduce the number 
of possible configurations available to the system, compared to the free boundary case, 
making it harder to excite the system, and tending to lock in the low temperature phase. 
This is observed in the phase diagrams where in the periodic boundary case the critical 
temperature estimates converge to the thermodynamic limit from above, whilst for the 
open boundary case the critical temperature estimates converge from below (see also 
table In table ^ we compare the critical-temperature estimates for the two boundary 
conditions for various key values of a and we find good agreement between the boundary 
conditions. 

The two estimated phase diagrams are, however, quite different at first sight, 
notably for a less than about —1.3. Before trying to reconcile the differences, let us 
summarise the similarities. For a = —1 there exists a finite-temperature phase transition 
with Tc = 0.37 ± 0.01. In table ^ we show finite-size estimates for the exponent u 
calculated using the two different boundary conditions. The extrapolated value of u 
for the open boundary conditions at a = — 1 indicates a transition in the three-state 
Potts universality class. The finite-size estimates for a = —1 with periodic boundary 
conditions are still too far from the thermodynamic limit to permit good extrapolation. 

The point a = — 1 with T = corresponds to a special point on the phase 
diagram, since it is the point which separates the ferromagnetic ground state [s = 0) 
from the frustrated highly degenerate ground state [s finite). It therefore corresponds 
to a transition point, which one would expect to be connected to the rest of the 
phase diagram. The phase diagram with open boundary conditions clearly shows a 
reentrant phase ending in a multiphase point which we identify with the special point 
T = 0, a = — 1. Whilst the periodic boundary results do show the begining of a 
reentrant phase, the results are less clear. However, as we have seen, there is good 
agreement between the extrapolated critical temperatures right up to « = —1. The 
transition lines become progressively more rounded as the system width increases. The 
failure of the transition lines to reach T = for the finite widths is probably related to 
a stabilisation of the low-temperature phase by the periodicity condition. To verify this 
conjecture, we compare the exact phase diagram found by Andre et al |^ with both 
periodic boundary condition and open boundary condition calculations, see figure 0. 
The convergence is excelent in both cases but, as for the Potts model, the periodic 
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boundary condition results fail reach T = (the transition lines for widths 12 and 
14 reach T = 0.444), whereas the open boundary condition transition lines do extend 
right down to T = 0. In the light of these considerations it seems to us likely that 
the point T = 0, a = —1 remains a multiphase point in the Potts model, where the 
paramagnetic, ferromagnetic and partially ordered phases meet. This view is conforted 
by the observation that the residual entropy per site at zero temperature is higher for 
a = — 1 than for either a < —1 or a > —1. This means that the ground state is 
more disordered at a = — 1 than for a < —1 presumably due the the coexistence of the 
paramagnetic and partially disordered phases. 

In models with competing interactions, as is the case here, one may define disorder 
lines which indicate where the nature of the correlation function changes from monotonic 
(due to the ferromagnetic interactions) to oscillatory (due to the antiferromagnetic 
interactions) fl^ . It is easy to show that the correlation function (in the transfer 
direction) is related to the first and second largest eigenvalues in absolute value (the 
largest being always positive) through: 



If Ai is positive, then G{x) is monotonically decreasing, but if Ai is negative then 



and is oscillatorily-decreasing. If we define A^'^ as the largest positive eigenvalue, and 
Xi as the largest (in absolute value) negative eigenvalue, then the disorder line is found 
when A^^ = |Aj"|. 

The disorder lines for various lattice widths are shown in figure ^ as lines. It 
is common to find disorder lines associated with reentrant phases. In this case the 
disorder lines must enter the cusp caused by the reentrant paramagnetic phase; since 
the disorder line corresponds to a change in behaviour of the correlation function it may 
not cross a critical line. The disorder lines for the two different boundary conditions 
may be seen to converge to a common line as the thermodynamic limit is approached. 
The limiting disorder line is curved, contrary to the Ising case, see [Appendix A| , were 
the disorder line is calculated analytically and found to be a vertical line at a = — 1. 
The two transition lines join in the form of a cusp, with the disorder line between them 
(as in the Ising case). The form of the disorder line may then be seen as confirming the 
existence of the reentrant paramagnetic phase. 

Phenomenological renormalisation indicates two possible phase diagrams depending 
on which boundary conditions were taken: Either the transition line extends to a — cxd 
and the phase diagram resembles qualitatively a distorted version of the Ising model 
Phase diagram, or there is a second reentrance. The true thermodynamic phase diagram 
may correspond to one of the two posibilities discussed, or to a combination of the two, 
or indeed to neither. 




(12) 
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In order to shed some light on this issue it is interesting to look at a variety of 
average thermodynamic quantities, such as the magnetisation, energy, average numbers 
of each type of bond that are satisfied, etc. It is possible to calculate with relative ease 
such local average quantities in the framework of a transfer matrix calculation. 

In order to calculate the magnetisation it is necessary to first calculate the 
probability of having configuration in column x. This probability is simply the 
ratio of the partition function restricted to having configuration and the unrestricted 
partition function, which in terms of transfer matrices may be written: 

Tr I C^)(C^ I T^-^ I 
PiC.) = ^ '^/^^ ^- (13) 

Writing \Cx) in terms of the eigenvectors of T and taking N to infinity and using the 
fact the T is a real symmetric matrix leads to: 

p{C.) = (0|C.)^ (14) 

The magnetisation for the Potts model is defined as: 

m = (g(nmax) - 1) , (15) 

q-1 

with nmax the maximum of {pi}, where pi are the densities of spins in each of the Potts 
states. Defining m{C) as the magnetisation for a column state as 

1 / iV^ax(C) .\ .... 

m{C) = Iq 11, (16) 

with A^max the largest of the numbers of each spin type present in configuration C. The 
magnetisation is then given as 

m = ^ m{C)p{C) 
c 

= Y^Hcmcf. (17) 

c 

Since (0|C)^ is the probability of finding a column in configuration C, then the 
probability of finding a column with magnetisation m is given by 

P{m) = Y.6{m{C) -m){0\Cf. (18) 
c 

It is important to note that this is not the probability of a magnetisation m for the 
system as a whole. Since this probability distribution is restricted to a column it gives 
a finer probe into the probable column configurations. This is plotted for a = —1.3 and 
T = 0.2 as a probability density function in figure §by dividing each probability by the 
interval between allowed values of the magnetisation for a given lattice width. It may 
be clearly seen that the probability density is developing two peaks, one at m = 1 and 
a second at a substantially lower value of m. This implies that in the low temperature 
phase each column has a high probability of being ordered ferromagnetically (along the 
y direction)or of being (partially) disordered. How these are arranged in the x direction 
is not easy to calculate in the transfer matrix formalism, but preliminary Monte Carlo 
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calculations ||18|| seem to indicate that the system alternates in the x direction between 
ferromagnetically ordered and disordered columns. The probability distribution shown 
is virtually identical for the two considered boundary conditions and for a = —1.3 and 
a = —2. The peak at m = 1 is lost in the high temperature phase, clearly indicating 
the absence of magnetic order. No evidence of two different frustrated phases could be 
found from the distribution of the magnetisation probability. 

To calculate quantities such as the energy per spin or the densities of satisfied bonds 
it is necessary to calculate the probabilities of having a configuration Cx in column x 
and Cx+i in column x + 1 since these quantities are defined in terms of pairs of spins, 
which may belong to two adjacent columns. The calculation of the joint probability 
p{Cx,Cx+i) follows exactly as above: 

piCx, = ^ Tr{T^} ^^^^ 

Defining nferro(Cx, C^+i) and ^iantifcrro(Cx7 Cx+i) as the densities of satisfied bonds, we may 
calculate the average values: 

(^ferro) ^ ^ '^ferro iCx,Cx+i)piCx,Cx+i)- (21) 

and similarly for (riantiferro)- 

These two quantities are plotted in figure |^ for (A) a = —1.3 and (B) a = —2 
with periodic boundary conditions. It may be seen that there is a clear difference of 
behaviour between a = —1.3 and a = —2. In the first case the two densities may 
be seen to drop away from their maximal values at the different temperatures (shown 
approximately by arrows in figure ^ (A)). In the second case this same drop off occurs at 
the same temperature for the two different densities. The results for the densities show 
a high degree of agreement between the two different boundary conditions. Both show 
a qualitative difference of behaviour between a = —1.3 and a = —2, and evidence of 
possible phase transitions for both values of a, indicated by arrows in figure |^, two for 
a = —1.3 and one for a = —2. The specific heat does not diverge with system size for 
either of these values of a, but the peaks of the specific heat converge to a temperature 
consistent with the critical temperature shown in figure ^ (A) for a less than about —1.5 
and all a > — 1. 

We conclude that there is sufficient evidence to conjecture a phase diagram as shown 
semi-schematically in figure [ly, with a possible second partially ordered phase at low 
temperature, (PO II). This additional phase, compared to the piled- up-domino Ising 
model, is only presented speculatively; the width of this phase reduces as system size 
increases (see figure |^) and may disappear in the thermodynamic limit. If such a phase 
were to exist then it is likely to terminate at a multicritical point (A). The presence of 
such a multicritical point could provide evidence of such a phase. 
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5. Conclusions 

The main results presented in this article relate to the Potts piled-up-domino model, 
where we have presented behaviour which is remarkably different from the equivalent 
Ising case, notably a finite-temperature transition and reentrant behaviour. It is unusual 
to see reentrance in such a simple model. Models showing reentrance normally contain 
competing nearest-neighbour and next-nearest-neighbour interactions |^ or interaction 
between different spin types [^. The additional entropy available due to the greater 



number of spin states stabilises the ferromagnetic ground-state as the temperature is 
increased (with a = —1), contrary to the naive anticipation that the extra degrees of 
freedom will disorder the system more. This may be understood as follows: in the 
Ising case the ground states may be mapped one to one onto dimer coverings of the 
lattice, but when the g-state Potts model is considered the mapping is no longer one 
to one, and the additional degeneracy is not uniformly distributed. This reweighting 
must apply to all the low lying states, and must serve to stabilise the ferromagnetic 
phase as the temperature is increased through an order by disorder mechanism. In the 
zig-zag case the different dimer coverings have the same weight in the ground state, as 
in the Ising case, and we do not observe an equivalent stabilisation of the ferromagnetic 
ground state. 

The exact nature of the phase diagram and the low temperature phases for a < — 1 
remain to be fully determined, and will doubtless require other methods. This will be 
the object of a future study. 
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Appendix A. Disorder line for the Ising piled-up-domino model 

The free energy and phase diagram for the Ising piled-up-domino model was solved 
exactly by Andre a/ 0, but to the best of our knowledge the disorder line has not 
been found for this model. 

Andre et al give the (dimensionless) free energy per spin as: 

/ = log2 + — dg^dgj/log A(g^, qy), (A.l) 



with 



^{<lx,qy) = ^[cosh2/3Jcosh2a/3J + cosh^2/3Jcosh2(l + a)/3J 

- sinh^ 2pJ cos q-y — 2cosh2/?Jsinh2(l + a)f3Jcosqx (A. 2) 
-|- sinh 2/5Jsinh 2l3aJ cos2qx], 
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from which we deduce the form of the generic correlation function as 
leading to 

The modulation of (^(x, ?/) in the x {y) direction is governed by the maximum of 
G{qx,qy), or the minimum of A{qx,qy), with respect to qx (qy). Since the modulation 
is anticipated in the x-direction (see[|l), we seek the turning points of A{qx,qy) with 
respect to q^: 

OA 

—— = 2cosh2/?Jsinh2(l+a)/?Jsing^-2sinh2/3Jsinh2a/3Jsin2ga.,(A.5) 
oqx 

which becomes zero when qx = 0, vr or when 

_ cosh2/3Jsinh2(l + a)/3J 

2sinh2/5Jsinh2a/3J ' ^ ' 

Checking the second derivative of A it may be seen that the q'^ = solution is a 
minimum as long as 

cosh2/?Jsinh2(l + a)/?J> 2 sinh2/?Jsinh2a/3J, (A.7) 

and the q^ = vr solution is a minimum as long as 

cosh2/5Jsinh2(l + a)/3J> -2 sinh 2/5Jsinh 2a/3J. (A.8) 

These two conditions correspond exactly to the boundaries within which the solution 
( |A.(j| ) exists, which must correspond to a maximum. 

In the range of values of a for which the two solutions qx = and = tt are 
minima it remains to determine which corresponds to the absolute minimum. This may 
be determined by considering 

A{7f, qy) - A{0, qy) = 4cosh2/5Jsinh2(l + a)/3J. (A.9) 

This changes sign when a = — 1. For a > —1 the leading mode of the correlation 
function is the q^ = mode, while for a < — 1 the leading mode is the = mode. 
This gives the location of the disorder line as a = — 1 for all temperatures. 
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Figure 1. The piled- up-doniino model (A) and the zig-zag model (B). The interaction 
energy is Ji along the solid lines and J2 along the dashed lines. 
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Figure 2. Construction of the transfer matrix by single spin addition 
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Figure 3. Plaquettes for the construction of the zig-zag model. Solid lines show 
interactions of strength Ji and dashed lines interactions of strengh J2. The arrow 
shows path to prefixed spin in the case that the marked bond is frustrated. 
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Figure 4. Phase diagram for the Potts Zig-Zag model 
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Figure 5. Finite size estimates for the exponent v for the Potts zig-zag model 
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Figure 6. Phase diagram for the Potts piled-up-domino model with: (A) periodic 

boundary conditions and (B) free boundary conditions. The disorder Hnes are shown 
as hnes, and the phase boundaries as symbols. The numerical uncertainty is smaller 
than the size of the symbols. 
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Figure 7. Exact phase diagram for the Ising Piled-up-domino model, found by Andre 
et al ||] shown as a solid line along with the finite size estimates calculated fi'om the 
phenomenological renormahsation group for widths 6 and 8 (O 8 and 10 (□), 10 and 
12 (<C>) and 12 and 14 (A). In (A) periodic boundary conditions are used, whilst in 
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Figure 9. The set of lines in each graph labeled (1) shows the number of satisfied 

antifcrromagnetic bonds normalised by the total number of antifcrromagnctic bonds, 
while the set of linos labeled (2) shows the number of satisfied ferromagnetic bonds 
per spin. Figure (A) corresponds to a = —1.3 and figure (B) to a = — 2 both with 
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Figure 10. Possible thermodynamic phase diagram for the three-state Potts piled- 
up-domino model. The disorder line is shown as a dashed line, and the dotted line is 
a possible extra transition line between two partially ordered phases (see text) with a 
possible multicritical point, marked as A on the phase diagram. 



Critical behaviour of fully-frustrated Potts models 
Tables and table captions 



Table 1. Ground-state entropy for the Potts pilcd-up-domino model 



r 


ni <r' 1 


n ^ 1 


— — 1 


n- — — 1 




Free b.c.'s 


Periodic b.c.'s 


Free b.c.'s 


Periodic b.c.'s 


2 


0.5025263 


0.5025263 


0.635098 


0.6350983 


4 


0.3854409 


0.3588588 


0.471894 


0.4333004 


6 


0.3521346 


0.3238563 


0.430075 


0.3934622 


8 


0.3369267 


0.3113484 


0.412780 


0.3825248 


10 


0.328266 


0.3057865 


0.403783 


0.3789744 


12 


0.3226777 


0.3028948 


0.398421 


0.3776883 


CX5 


0.2999 


0.2972 


0.3814 


0.3766 


Ising 


ilog((l + V5)/2) 




G 




Exact 


= 0.240659- •• 




0.2916- • • 





Table 2. Critical temperature estimates for the zig-zag model. 



L/L' 


a = 1 


a = 


4/6 


1.0028517 


0.679435201 


6/8 


0.997615939 


0.675820268 


8/10 


0.996228353 


0.674737229 


10/12 


0.99569457 


0.674343044 


oo 


0.995 


0.674 


Exact 


l/log(l -1-73)0 


(1) 




= .994973 • • • 


= .6737596022 • • • 



(1) Critical temperature for a = solution of e"^/^ — 3e^^^ — 6e^/^ — 1 = 0. 
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Table 3. Critical temperature estimates for the piled-up-domino model 



r It' 

Li 1 Lj 


a = 


: 




a = -1 




a = 


-1.3 


PBC 


FBC 


PBC 


PBC-y 


FBC 


PBC 


FBC 


2/4 


0.764578 


0.623342 


0.482104 


0.428436 






0.081726 


4/6 


0.739625 


0.677582 


0.456832 


0.416958 


0.272384 


0.350795 


0.189215 


6/8 


0.727985 


0.696787 


0.417748 


0.393249 


0.310488 


0.329505 


0.235781 


8/10 


0.724292 


0.705696 


0.395310 


0.381282 


0.328203 


0.315065 


0.253698 


10/12 


0.722934 


0.710524 


0.383267 


0.367335 


0.338479 


0.305799 


0.262133 


oo 


0.722 


0.722 


0.361 


0.367 


0.379 


0.283 


0.274 


exact 


l/log(4) = . 


.721358- •• 













PBC: Periodic boundary conditions 
FBC: Free boundary conditions. 

PBC-y: Periodic boundary conditions with transfer direction perpendicular to dashed 
lines (see Figure |l|) 



Table 4. Estimates for the exponent v for the piled-up-domino model 



Periodic Boundary Conditions Free Boundary Conditions 
L/L' 





a = l 


a = 


a = —1 


a = l 


a = 


a = -1 


2/4 


0.807155 


0.727115 


0.457417 


0.897734 


0.770501 




4/6 


0.827697 


0.754573 


0.505535 


0.891758 


0.836442 


2.049413 


6/8 


0.837053 


0.792602 


0.573436 


0.883714 


0.852665 


1.354167 


8/10 


0.839283 


0.811782 


0.639087 


0.877283 


0.857856 


1.162991 


10/12 


0.840722 


0.821423 


0.697795 


0.872277 


0.859395 


1.075589 


CO 


0.842 


0.837 


1.139 


0.836 


0.860 


0.872 



